Local moments and magnetic order in the two-dimensional Anderson-Mott transition 
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We study the role of electronic correlation in a disordered two-dimensional model by using a 
variational wave function that can interpolate between Anderson and Mott insulators. Within this 
approach, the Anderson-Mott transition can be described both in the paramagnetic and in the 
magnetic sectors. In the latter case, we find evidence for the formation of local magnetic moments 
that order before the Mott tran sition. Th e charge gap opening in the Mott insulator is accompanied 
by the vanishing of the lim 9 ^o { n q){'n-q) (the bar denoting the impurity average), which is related 
to the compressibility fluctuations. The role of a frustrating (second-neighbor) hopping is also 
discussed, with a particular emphasis to the formation of metastable spin-glass states. 

PACS numbers: 



The combined action of electron-electron interaction 
and disorder is known to heavily influence the physi- 
cal behavior of electron systems. [l| Recently, the ob- 
servation of metallic behavior in high-mobility two- 
dimensional electron-gas devices [2j has opened new per- 
spectives in this subject, suggesting the possibility that a 
metallic behavior could be stabilized by a strong electron- 
electron interaction in two dimensions, in spite of the 
standard scaling theory of Anderson localization. 0, Q 
Such a proposal was first put forward theoretically by 
means of a weak-coupling renormalization group ap- 
proach within a Fermi-liquid description, [f| and later 
developed along similar directions. [6|, [7j A common fea- 
ture of the above renormalization-group calculations is 
the crucial role played by the spin fluctuations that grow 
large as the renormalization group procedure is iterated. 
This tendency, which has been interpreted as signaling 
the emergence of local moments, suggests that electron- 
electron correlations become effectively very strong that, 
in turn, makes doubtful the validity of a Fermi liquid 
description. [H, Q 

Apart from the debated issue of a metal-insulator tran- 
sition in two-dimensional high-mobility devices, HI 
Hil . Fuj i there are less controversial systems where the 
role of strong correlations concomitantly with disorder 
is well testified. Particularly emblematic is the case of 
Si:P and Si:B, 3, IH which are three-dimensional ma- 
terials that show a bona fide metal-insulator transition. 
Here, the randomly distributed impurities form a very 
narrow band within the semiconducting gap. Since the 
local Coulomb repulsion is sizable compared to the width 
of the impurity band, this system is particularly suit- 
able to investigate the interplay between disorder and 
interaction. Indeed, clear signatures of local magnetic 
moments are found in several thermodynamic quanti- 
ties. 0, 17, [ID, [l9[ Theoretically, the interplay of dis- 
order and interaction is a very difficult question. Any 
approach based on a sin gle-p article description, like un- 
restricted Hartree-Fock, [2fJ, [21| can uncover the emer- 



gence of local moments only if spin-rotational symmetry 
is explicitly broken, introducing spurious effects due to 
magnetism that can be dealt with using further approx- 
imate schemes. @, [24| More sophisticated approaches, 
like those based on dynamical mean-field theory, 1251 1 can 




in principle manage without magnetism, 
they usually miss important spatial correlations. 

In this Letter, we will generalize the variational ap- 
proach that has been successfully used to describe the 
Mott transition in finite-dimensional clean systems [29l . 
I30I l3lj ]. We will show that, for a half-filled disordered 
Hubbard model on a square lattice and when the vari- 
ational wave function is forced to be paramagnetic, the 
Anderson to Mott insulator transition exists and it is con- 
tinuous. When magnetism is allowed, we find two succes- 
sive second order phase transitions: from a compressible 
paramagnetic Anderson insulator with local moments to 
a compressible magnetic Anderson insulator and then to 
an incompressible magnetic Mott insulator. Unlike pre- 
vious unrestricted Hartree-Fock [2lJ or Monte Carlo cal- 
culations, [22| we do not find any evidence of an inter- 
mediate truly metallic behavior. 

We consider a half-filled Hubbard model on a square 
lattice with on-site disorder: 



+ H.c. + '^2(e i m + U ni^n^i), (1) 



where c\ a (ci j(T ) creates (destroys) one electron at site 

% with spin er, n il<7 = c\ a c i<y , and = ^ V' e i 
are random on-site energies chosen independently at each 
site and uniformly distributed in [—D,D]. tij are the 
hopping parameters that we will consider limited either 
to nearest-, tij — —t, or to next-nearest-neighbor, tij = 
—t', sites. In the calculations we will consider 45 degree 
rotated clusters with N = 2n 2 sites, n being an odd 
integer, and periodic boundary conditions, so that the 
non-interacting ground state is always non-degenerate at 
half filling. 



Following the approach developed for clean sys- 
tems, p9| we define a variational wave function contain- 
ing a Gutzwiller and a long-range Jastrow factor that 
apply to an uncorrelated state: 



q,co = 



m = v G j\4> ), 



(2) 



where |$o) is the ground state of a non- interacting Hamil- 
tonian with the same hopping parameters as in Eq. ([1]) 
but with variational spin-dependent on-site energies 1i a 
to be determined by minimizing the total energy. A para- 
magnetic wave function is obtained by forcing e^j = e^, 
while, to discuss magnetism, we allow the wave func- 
tion to break spin-rotational symmetry with e< + ^ e< 1. 
T'g — ex P 3i n i ] i s a Gutzwiller correlator that de- 
pends upon the site- dependent parameters <7j's, while 



J 



exp 



1 



is a Jastrow 



factor. The latter one spatially correlates valence fluctua- 
tions, Sni = (n, — 1} 7^ 0, on different sites, binding those 
with 5ni Snj < and unbinding those with Sni Snj > 0. 
This fact has been shown to be crucial to describe a Mott 
transition in clean systems. 2^, 3(3] We shall assume that 



Vij is translationally invariant, which makes the numer- 
ical calculations feasible but neglects any clustering ef- 
fects. All the parameters contained in the variational 
wave function j^), i.e., £i jCr , gi, and v^j, [13] are opti- 
mized to minimize the variational energy by using the 
Monte Carlo technique of Re f. 1331 ] . 

As discussed in Refs. [2^, (3C| for clean systems, it is 
possible to discriminate variationally metals from Mott 
insulators by looking to the equal-time density-density 



structure factor N„ 



I where n q is 



the Fourier transform of the electron density rij. Indeed, 
N q ~ \q\ implies the existence of gapless modes, while 
N q ~ \q\ 2 indicates that charge excitations are gapped. 
Moreover, there is a tight connection between the long- 
wave-length behavior of N q and the Fourier transform 
of the Jastrow factor v q , namely v a ~ 1 / 1 <7 1 for a metal 
and v q ~ 1/|?| 2 for an insulator. [29|, |30( This distinc- 
tion should equally work in ([1]) after disorder average. 
However, particular care must be taken to interpret N q 
in a disordered system, where the structure factor in- 
cludes a disconnected term, N q lsc = (n q )(n- q ) (where 
the quantum average is taken at fixed disorder config- 
uration and the overbar indicates the disorder average) 



as well as a connected one, i.e. 



7yconn 
9 



N„ 



l\ q 



For a clean system, the disconnected term gives rise to 
the elastic scattering peaks at q equal to the reciprocal 
lattice vectors, the Bragg reflections. On the contrary, 
in the presence of disorder iV^ lsc is finite for any finite 
momentum q. 34] The diagrammatic representation of 
N q lsc is shown in Fig. Q] and one can realize that, for 
q —> 0, it reduces to the electron compressibility fluctu- 
ations. For non- interacting electrons, iV^ lsc is finite for 
q — > 0, whereas N q onn ~ \q\, indicating the absence of a 
gap in the spectrum of charge-density fluctuations. 34 ] 




FIG. 



1: (Color online) Density-density fluctuations N q 1BC = 
(n q ){n- q ) . Dotted lines denote impurity averages, and the 
squares indicate vertex corrections that include both inter- 
action and impurity insertions. Continuous lines are fully 
corrected Green's functions. 
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FIG. 2: (Color online) (a) Connected term of the density- 
density correlation function N q onn divided by |qr|. (b) Jastrow 
parameters v q multiplied by \q\ 2 . (c) Disconnected term of 
the density-density correlation function iV^ lsc as a function 
of U. (d) Fluctuations of the on-site variational energies A e 
and of the local densities. All calculations have been done for 
D/t = 5. 



We start our analysis with the case of nearest-neighbor 
hopping only by using a paramagnetic wave function, 
namely imposing e^j = e^.j. In Fig. [21 we show the varia- 
tional N q onn and the Fourier transform of the optimized 
Jastrow potential v q for different values of the interac- 
tion U and D/t = 5 (we take such a large value of D 
in order to have a localization length that, at U = 0, is 
smaller than the numerically accessible system sizes). A 
clear change in the behavior of these quantities is ob- 
served at U™ ott /t = 11.5 ± 0.5. For small values of 
the electron interaction, N q onn 



q\ and v q ~ l/\q\, 



whereas N^ onn 



\q\ and v q ~ VM * n ^ ne strong- 



coupling regime. The latter behavior is symptomatic of 
the presence of a charge gap hence of a Mott insulating 
behavior. J2i| We notice that, for the clean case D = 
and within the same approach, a metal-insulator transi- 
tion at c/ c Mott = 8.5 ± 0.5 was found, [HI indicating that 
disorder competes with U and pushes the Mott transi- 
tion to higher values of U/t. It should be emphasized 
that, with respect to the clean system, for U < [/Mott ^ 
Nq° nn ~ \q\ is not associated to a metallic behavior but 
only to a gapless spectrum, also characteristic of an An- 



3 



derson insulator. Remarkably, we find that the Mott and 
Anderson insulators can also be discriminated through 
the behavior of the lim g ^ N^ lsc . In Fig. [5] we plot this 
quantity for different values of U , demonstrating that it 
is finite in the Anderson insulator, whereas it vanishes in 
the Mott phase. This identifies a simple and variation- 
ally accessible order parameter for the Anderson-Mott 
transition. 

Even though within this approach we cannot access 
dynamical quantities like DC conductivity, hence we 
can not address the question of a possible stabilization 
of a conducting phase with moderate Coulomb repul- 
sion, [22] we note that the linear slope of N^ onn has a 
non-monotonic behavior as a function of U, showing a 
peak for U/t ~ 7 that indicates an accumulation of low 
energy states around the Fermi energy. The same quali- 
tative behavior is also present in the fluctuations of the 
local densities, 5n 2 = l/JV^.i^-) - K) 2 ). Though the 
single-particle eigenstates of the variational Hamiltonian 
may have a very long localization length, because of the 
suppression of the effective on site disorder ej, yet this 
length is still finite in two dimensions hence the many- 
body wave function | 1 E r ) always describes an Anderson 
insulator below the Mott transition. Indeed, as shown in 
Fig. [2 the fluctuations of the on-site variational disorder 
A 2 = 1/N J2i e? J2i e ~i) 2 are always finite, though 

sizably renormalized by the electron interaction U. 

Let us now move to the more interesting case in which 
we allow magnetism in the variational wave function, 
which amounts to permit e^j ^ e^j. In this case the 
ground state may acquire a finite local magnetization 
on each site m, = — Tii t \. A magnetically ordered 
phase will have a finite value of the total magnetization 
M = 1/N^2j e lRjt ^m,j for a suitable momentum Q, like 
for instance Q = (tt, tt) for the Neel state. In the pres- 
ence of disorder, a finite value U^ F is needed to have 
long-range antiferromagnetic order. We find that, also 
in presence of a small t', U^ F < C/^ ott , giving rise to 
an extended region with antiferromagnetic order and fi- 
nite compressibility (i.e., a vanishing charge gap). These 
results are i n ag reement with previous mean-field calcu- 
lations. [U l36l . [37l | In Fig. [3J we show the results for 
t' = cither by fixing D/t = 5 and varying U (for which 
Uf/t = 6.5 ± 0.5 and t/ c Mott /i = 10.5 ± 0.5) or by fixing 
U/t = 4 and changing D (for which Df ott /t = 1 ± 0.5 
and D^ F jt = 2.5 ± 0.5). We note that the onset of anti- 
ferromagnetism is preceded by a magnetically disordered 
phase (i.e., M = 0) in which local moments appear. In 
Fig. [H the pattern of the local density (rn) and local 
magnetization (mi) are shown for a typical realization of 
disorder. For U/t = 4, the ground state is an Ander- 
son insulator with a large number of empty and doubly 
occupied sites with m, ~ 0. However, some sites have 
finite magnetization, but they are not spatially corre- 
lated hence long-range magnetism is absent. We interpret 
these magnetic sites as local moments. When the electron 
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FIG. 3: (Color online) Staggered magnetization M for Q = 
(tt,tt) and compressibility fluctuations N q ' BC as a function of 
U for disorder D/t = 5 (upper panel) and as a function of D 
for U/t = 4 (bottom panel). Fluctuations of the on-site vari- 
ational energies A,; and of the local densities (middle panel) . 
Calculations have been done for N — 98 and error-bars indi- 
cate the average over different realizations of disorder. 
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FIG. 4: (Color online) Local density (m) (upper panels) and 
local magnetization (mi) (lower panels) for a given disorder 
realization with D/t = 5 and different values of U/t. The 
black contour shows the elementary cell of the lattice which 
it is repeated to mimic the infinite lattice with periodic bound- 
ary conditions. 



interaction U increases, the number of magnetic sites in- 
creases rapidly and the local moments eventually display 
the typical staggered pattern of Neel order. Nevertheless, 
charge excitations are still gapless, with N^ onn ~ \q\. For 
U/t = 12 the system is a gapped insulator with antifer- 
romagnetic order and a vanishing compressibility. Vari- 
ationally, the charge gap opens by the combined effect 
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FIG. 5: (Color online) Local magnetization (m.) for the best 
variational state (lower panels) and for a metastable solution 
(upper panels) for a given disorder configuration with D/t = 5 
and t'/t = 1. 



of the Jastrow correlations, i.e., v q ~ l/M 2 , and the fi- 
nite antiferromagnetic gap in the mean-field Hamiltonian 
(due to staggered e^'s). 

In the presence of a large frustrating hopping t'/t > 0.9 
we find evidences of a spin glass behavior. In the large 
U regime, the optimal wave function displays magnetic 
long-range order with Q = (ir, 0) or (0, n). However, the 
energy landscape contains other local minima very close 
in energy in which most of the sites of the lattice have a 
net magnetization but an overall vanishing magnetic or- 
der, a "glassy" spin patterns, see Fig. [5] These solutions 
are incompressible, i.e., iV^ lsc ~ and, therefore, may 
be viewed as disordered Mott insulators. By decreasing 
the interaction strength U, these metastable states turn 
compressible, still having a large number of local mo- 
ments. However, the actual variational minimum shows, 
as before, a Mott transition from a Mott to an Anderson 
insulator, both magnetically ordered, followed, at lower 
U , by a further transition into a paramagnetic Anderson 
insulator. The only role of t' is to shrink the region in 
which a magnetic Anderson insulator is stable. 



In conclusion, we have shown that a relatively sim- 
ple variational wave function is able to describe the 
Anderson-insulator to Mott-insulator transition in two 
dimensions. In the paramagnetic sector, this phase tran- 
sition is continuous, in agreement with dynamical mean 
field theory. [2f 



27] When spontaneous spin symmetry 



breaking is allowed, we find two successive phase tran- 
sition, the first from a paramagnetic Anderson insula- 
tor to a magnetic one, followed by a transition from a 
magnetic Anderson insulator to a magnetic Mott insu- 
lator. Upon increasing frustration, the stability region 
of the magnetic Anderson insulator decreases. In gen- 
eral, the paramagnetic Anderson insulator develops local 
magnetic moments, but we do not find any evidence of a 
truly metallic behavior induced by interaction. 
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